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Abstract 

The production of beauty quarks in ep interactions has been studied with the 
ZEUS detector at HERA for exchanged four-momentum squared > 10 GeV^, 
using an integrated luminosity of 363 pb~^. The beauty events were identified 
using electrons from semileptonic b decays with a transverse momentum 0.9 < 
< 8 GeV and pseudorapidity |?7^| < 1.5. Cross sections for beauty production 
were measured and compared with next-to-leading-order QCD calculations. The 
beauty contribution to the proton structure function F2 was extracted from the 
double-differential cross section as a function of Bjorken-x and Q^. 
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1 Introduction 



The production of heavy quarks in ep colhsions at HERA is an important testing ground 
for perturbative Quantum Chromodynamics (pQCD), since the large 6-quark mass pro- 
vides a hard scale that allows perturbative calculations to be made [I1I2]- The dominant 
production process is boson-gluon fusion (BGF) between the incoming virtual photon and 
a gluon in the proton. Beauty production has been measured using several methods by 
the ZEUS P4TT] and the HI [T214T8] collaborations both in deep inelastic scattering (DIS), 
i.e. for large exchanged four-momentum squared, Q^, and also in photoproduction, i.e. for 
~ GeV^. The measurements are reasonably well described by next-to-leading-order 
(NLO) QCD predictions. 

Most of the previous measurements of 6-quark production used muons to tag semileptonic 
decays of the B hadrons. This paper reports a measurement of beauty production in DIS 
using the semileptonic decays to electrons, 

ep-^e'bbX ^e'eX', 

in the kinematic range > lOGeV^. Using the electron channel allows a measure- 
ment of the decay leptons at lower transverse momentum and provides a complementary 
measurement, with independent systematics. 

An analysis of the same process in the photoproduction regime, based on data taken in 
1996-2000 (120 pb~^), used a likelihood-ratio test to extract the signal of beauty and 
charm semileptonic decays to electrons A similar method, adapted to the different 
kinematics of the DIS regime, was used for the measurement reported here. The analysis 
also benefited from improved tracking in the more recent data, which allowed the measured 
decay length of weakly decaying B hadrons to be used. 

In this analysis, the total visible cross section, a^^e, and differential cross sections as a 
function of Q^, the Bjorken scaling variable, x, the transverse momentum, p^, and the 
pseudorapidity of the electron, 77^, were measured. They are compared to a leading-order 
(LO) plus parton-shower (PS) Monte Carlo prediction and to an NLO QCD calcula- 
tion. The beauty contribution to the proton structure function F2, denoted as F^^, was 
extracted from the double-differential cross section as a function of and x and is 
compared with theoretical calculations. 

2 Experimental set-up 

This analysis was performed with data taken from 2004 to 2007, when HERA collided 
electrons or positrons with energy = 27.5 GeV with protons of an energy of 920 GeV, 
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corresponding to a centre-of-mass energy y/s = 318 GeV. This data-taking period is 
denoted as HERA II. The corresponding integrated luminosity is (363 ± 7) pb~^. 

A detailed description of the ZEUS detector can be found elsewhere [12] . A brief outline 
of the components that are most relevant for this analysis is given below. 

In the kinematic range of the analysis, charged particles were tracked in the central track- 
ing detector (CTD) [20] and the microvertex detector (MVD) [21]. These components 
operated in a magnetic field of 1.43 T provided by a thin superconducting solenoid. The 
CTD consisted of 72 cylindrical drift chamber layers, organised in nine superlayers cover- 
ing the polar-angl^ region 15° < 6 < 164°. The MVD silicon tracker consisted of a barrel 
(BMVD) and a forward (FMVD) section. The BMVD provided polar-angle coverage 
for tracks with three measurements from 30° to 150°. The FMVD extended the polar- 
angle coverage in the forward region to 7°. After alignment, the single-hit resolution of 
the BMVD was 24 fim and the average impact parameter resolution of the CTD-BMVD 
system for high-momentum tracks was 100 /xm. 

To estimate the ionisation energy loss per unit length, dE/dx, of particles in the CTD |22] . 
the truncated mean of the anode-wire pulse heights was calculated, which removes the 
lowest 10% and at least the highest 30% depending on the number of saturated hits. 
The measured dE/dx values were corrected by normalising to the average dE/dx for 
tracks around the region of minimum ionisation for pions with momentum, p, satisfying 
0.3 <p< 0.4 GeV [23]. 

The high-resolution uranium-scintillator calorimeter (CAL) [21] consisted of three parts: 
the forward (FCAL), the barrel (BCAL) and the rear (RCAL) calorimeters. Each part 
was subdivided transversely into towers and longitudinally into one electromagnetic sec- 
tion and either one (in RCAL) or two (in BCAL and FCAL) hadronic sections. The 
smallest subdivision of the calorimeter is called a cell. The CAL energy resolutions, 
as measured under test-beam conditions, are a{E)/E = 0.18/ VE for electrons and 
a{E)/E = 0.35//E for hadrons, with E in GeV. 

The luminosity was measured using the Bethe-Heitler reaction ep — )■ ejp by a luminosity 
detector which consisted of independent lead-scintillator calorimeter [2S] and magnetic 
spectrometer [26] systems. 



The ZEUS coordinate system is a right-handed Cartesian system, with the Z axis pointing in the 
proton beam direction, referred to as the "forward direction", and the X axis pointing left towards 
the centre of HERA. The coordinate origin is at the nominal interaction point. The pseudorapidity 
is defined as 77 = — In (tan |), where the polar angle, 6, is measured with respect to the proton beam 
direction. The azimuthal angle, 0, is measured with respect to the X axis. 
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3 Monte Carlo simulation 



To evaluate the detector acceptance and to provide the signal and background distri- 
butions for the likelihood-ratio test, Monte Carlo (MC) samples of beauty, charm and 
light-flavour events were generated, corresponding to eighteen, two and one times the 
integrated luminosity of the data, respectively. The Rapgap 3.00 Monte Carlo pro- 
gram [27] was used to generate the beauty and charm samples. The CTEQ5L [28] parton 
density functions were used and the heavy-quark masses were set to = 4.75 GeV and 
rric = 1.5 GeV. To simulate radiative corrections, the events were passed through the 
Heracles 4.6 [22] program. An inclusive MC sample containing all flavours was gen- 
erated using Djangoh 1.6 [30] interfaced to Ariadne 4.12 [31], where the quarks were 
taken to be massless. The CTEQ5D [28] parton density functions were used. 

For the acceptance determination, the distribution in the signal MC was reweighted 
in order to correct for observed differences between the measured and simulated distribu- 
tions. The corrections varied from -|-10% at low to —30% at high Q^. The 5-hadron 
lifetimes were corrected for differences between the simulated values and the world-average 
values [52] . 

Fragmentation and particle decays were simulated using the Jetset/Pythia model [35] . 
The lepton energy spectrum from charm decays was reweighted to agree with CLEO 
data [M]. The generated events were passed through a full simulation of the ZEUS 
detector based on Geant 3.21 [35]. The final MC events had to fulfil the same trigger 
requirements and pass the same reconstruction program as the data. 

4 Theoretical predictions and uncertainties 

Next-to-leading-order QCD predictions were obtained from the HVQDIS [3H] program in 
the fixed-fiavour-number scheme (FFNS) [37]. More details about the calculation can be 
found elsewhere [1]. 

The 6-quark mass (pole mass) was set to nib = 4.75 GeV. The renormalisation and fac- 
torisation scales, /iij and fip, were chosen to be equal and set to /ir = /ii;' = 
The parton density functions were obtained from the FFNS variant of the ZEUS-S fit [38] 
using the same 6-quark mass as in the HVQDIS calculation. The value of as{Mz) was set 
to 0.105. 

The Peterson fragmentation function [32], with = 0.0035 [10], was used to produce 
beauty hadrons from the heavy quarks. The semileptonic decay spectrum was taken from 
the Pythia Monte Carlo. The contributions from prompt and from cascade decays. 
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b — )■ c(c) — e, including b t ^ e and b ^ J/ip ^ e^e , were taken into account in the 
effective branching fraction, which was set to 0.217 [52] . 

To estimate the uncertainty on the theoretical predictions, the 6-quark mass was varied 
in the range nib = 4.5, 5.0 GeV, and the scales fiR, fip were varied independently by a 
factor of two up and down. The parameter ef, was varied by ±0.002. The parton density 
functions were varied within the total uncertainties of the fit. The uncertainty on the 
NLO QCD prediction for the total cross section is +15% and —16%, where the dominant 
contribution originates from the variation of the mass and the scales. 

The HVQDIS calculations were also used to extrapolate the visible cross sections to ^2^''. 

5 Data selection 

Events were selected online with a three- level trigger [121 HI] using a combination of 
triggers, which required a scattered electron to be detected in the CAL and/or the presence 
of an electron candidate from a semileptonic decay. Further details on the trigger chain 
can be found elsewhere [12] . Offline, the reconstructed scattered electron was required to 
have an energy Eg' > 10 GeV. The Z position of the primary vertex had to be within 
|Zvtx| < 30 cm. 

The final state of the electron-proton collision, including the scattered electron, was re- 
constructed from energy-fiow objects (EFOs) [13] which combine the information from 
calorimetry and tracking, corrected for the energy loss in the detector material. Each 
EFO was assigned a reconstructed four-momentum, = (p^,py,p^, E'*). Jets were re- 
constructed from EFOs using the algorithm [H] in the longitudinally invariant mode 
with the massive recombination scheme [15]. 

The following cuts were applied to select DIS events: 

• the photon virtuality, Q^, must be above lOGeV^, where this variable and Bjorken-x 
were reconstructed using the double-angle method [16]; 

• 0.05 < y < 0.7, where the inelasticity, y, was reconstructed using the Jacquet-Blondel 
method [17] for the lower cut and the electron method jlH] for the higher cut; 

• 40 < {E — pz)tot < 65 GeV, reconstructed using the four-momentum of the final state; 
this selects fully contained neutral-current electron-proton events for which E — pz = 
2 -Eg = 55 GeV; 

• Pt/ Et < 0.7, where Pt and E^ are the transverse momentum and the scalar transverse 
energy of the final state. This cut was applied to reduce the charged- current and non- 
ep backgrounds. 
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In order to estimate the decay length of the B hadron, a secondary vertex was fitted using 
all good tracks assigned to the jet |1H]. Good tracks were defined by a minimal transverse 
momentum, px > 0.5 GeV, at least four hits in the MVD and three or more superlayers 
passed in the CTD. Vertices with x^/dof < 6 and a distance from the interaction point 
within ±1 cm in the X-Y plane and ±30 cm in the Z direction were taken. 

The decay length, d, was defined as the distance in XY between the secondary vertex 
and the interaction poinlj^, projected onto the jet axis. The sign of the decay length 
was assigned using the axis of the jet to which the vertex was associated; if the decay 
length vector was in the same hemisphere as the jet axis, a positive sign was assigned to 
it, otherwise the sign of the decay length was negative. Negative decay lengths, which 
originate from secondary vertices reconstructed on the wrong side of the interaction point 
with respect to the direction of the associated jets, are unphysical and caused by detector 
resolution effects. A small correction [12] to the MC decay-length distribution was applied 
in order to reproduce the data with negative values of decay length; 5% of the tracks in 
the central region were smeared and an additional smearing to tracks in the tails of the 
decay-length distribution was applied. 

Electron candidates from semileptonic decays of b quarks were selected from the EFOs 
having a transverse momentum, p^, satisfying 0.9 < < 8 GeV in the pseudorapidity 
range |?7'^| < 1.5, and consisting of a track matched to a single calorimetric cluster. To 
reduce the hadronic background, at least 95% of the EFO energy had to be deposited 
in the electromagnetic part of the calorimeter. Candidates in the angular regions corre- 
sponding to the gaps between FCAL and BCAL as well as between RCAL and BCAL 
were removed. To account for differences in the rj^ distribution in data and MC, the 
electron reconstruction efficiency in MC was corrected by 0.95 in the FCAL and RCAL 
regions and by 1.05 in the BCAL region. Electrons from identified photon conversions 
were rejected [TT] . 

The electron candidate was required to be associated with a jet using the following crite- 
ria: 

• the jet was required to have a reconstructed vertex of good quality as defined above; 

• the jet had to have p^^* > 2.5 GeV and |?7-''^*| < 2.0; 

• the distance AR = ^(r^jct _ ^e)2 _|_ (^^^ct _ 0e)2 < i q. 

• if there is more than one candidate jet, the jet closest in AR to the electron candidate 
was chosen. 

^ In the X-Y plane, the interaction point is defined as the centre of the beam ellipse, determined using 
the average primary vertex position for groups of a few thousand events, taking into account the 
difference in angle between the beam direction and the Z direction. The Z coordinate is taken as the 
Z position of the primary vertex of the event. 
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The combination of the momentum cut and the jet association reduces substantially the 
background from scattered electrons not identified as such. 

The main variable for the electron identification was dE/dx [7]. To reduce the major 
background of fake electrons in the candidate selection, a preselection cut was apphed on 
a likelihood-ratio test function j'^^/'^^ This function was calculated using dE/dx as 
discriminating input variable and testing the electron hypothesis. The distribution of this 
test function, as obtained from MC, for the particle types e^,7r^, p/p and is shown 
in Fig. [H The vertical line at —2\nTe^^'^^ = 3 indicates the cut, which rejects a large 
fraction of the background particles. 



6 Identification of electrons from semileptonic decays 

The electron candidates in the MC samples were classified into three different categories. 
The first category [b — )■ e) contains electrons from beauty decays, including direct semilep- 
tonic decays, cascade decays b — j- c(c) — )■ e, 6 — r — e and b ^ J /ip e+e". The second 
category (other e) contains all true electrons, which are not included in the beauty signal. 
These are mainly electrons originating from photon conversions, Dalitz decays, electrons 
from direct charm decays, or remaining DIS electrons. The third category (non-e) includes 
all candidates which are fake electrons. After the selection, the dominant contribution to 
the latter comes from pions, while the number of kaons or protons mimicking electrons is 
rather small. 

For the electron identification, the following three variables [7] were used as discrimi- 
nants: 

• dE/dx, as measured in the CTD; 

energy of the EFO as measured in the calorimeter, divided by the 
track momentum; 

• dccWi the depth of the central energy deposit within the CAL. 

The following discriminating variables were used to distinguish the origin of electron 
candidates: 

• the transverse-momentum component of the electron candidate relative to the 
direction of the jet axis. The shapes of the light-quark p^^^ distributions in the MC 
were corrected |12] using a background-enriched data sample. This variable is sensitive 
to b decays since electrons from b decays tend to have large p^^' due to the large b mass; 

• A0, the difference of azimuthal angles of the electron candidate and the missing trans- 
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verse momentum vector, defined as 




wliere is the negative vector sum of the EFO momentum transverse to the beam 
axis, 

and the sum runs over all EFOs. The variable A0 is sensitive to semileptonic decays 
of b and c hadrons due to the presence of the neutrino; 

• d/6d, the signed decay-length significance, where 6d is the uncertainty on d |^[50]. 
This variable is sensitive to the decay of c and b hadrons due to their long lifetimes. 

In contrast to the results of a previous ZEUS study [7], the separation power of A0 and 
p^r^^ is worse due to the lower jet momenta used here. Therefore it was not possible to 
separate the charm signal from the other particles in the electron background. 

Following the procedure of the previous study [7], the six variables were combined into 
one discriminating test-function variable, which is a ratio of likelihoods. For a given 
hypothesis of particle, i, and source j, the likelihood, Cij, is given by 

I 

where Vij{di) is the probability to observe particle i from source j with value di of a 
discriminant variable. The particle hypotheses i G {e,7i,K,p} and the sources, j G {6 
e, other e, non-e}, were considered. For the likelihood ratio test, the test function Tij was 
defined as 

= ' J '■J 
k,l 

The Oj, a'j denote the prior probabilities taken from MC. In the sum, /c, / run over all 
particle types and sources defined above. In the following, T is always taken to be the 
likelihood ratio for an electron originating from a semileptonic 6-quark decay, T = T^^b^e, 
unless otherwise stated. 



7 Signal extraction 

The combined MC sample was split into the three contributions as defined in the previous 
section. The beauty test function, T, was calculated separately for these three samples 
and for the data. The relative contributions of the three sources in the data, f^^^, 
fothl/ey fmm^e^y Were obtained from a three-component maximum- likelihood fit to the 
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T distributions. The fit range of tlie test function was restricted to —2 InT < 10 to remove 
the region dominated by background and where the test function falls rapidly. The for 
the fit is xVndf =18/28. 

The result of the fit is shown in Fig. [2] and corresponds to a scaling of the cross section 
predicted by the beauty MC by a factor of 1.32 ± 0.11. For the other two samples the 
scaling factors were determined to be ~1.1 for the electron background and ~1.3 for the 
non-e background. These factors were applied to the contributions shown in Figs. [3] and HI 

Figure [3] shows a comparison of the MC simulation to the data for the main variables 
used for the event selection. The Monte Carlo describes the data well. Figure H] shows 
the distributions for the variables in the likelihood-ratio test function, which are sensitive 
to the different origin of the electron candidates. In Figs. H] (a), (c) and (e), the three 
variables are shown for the selection used in the fit. Figures H] (b), (d) and (f), show the 
same distributions for a signal-enriched region, which is defined by a harder cut on the 
test function at —2 InT < 1.5. All distributions are reasonably well described. 



8 Cross-section determination 

The differential beauty cross section for a variable, was determined separately for each 
bin, k, from the relative fractions in the data obtained from the fit and the acceptance 
correction, Al"^^, calculated using MC events, 

dvfc Alte ■ ■ Avk 

where A^^ata ^^^q number of electron candidates found in the data bin, C is the inte- 
grated luminosity, Avk is the bin width and Cr is the QED radiative-correction factor. 
The acceptance is defined as 



rec 



A. -'^^ 



ATtrue ' 
b->e 



where N^^^ is the number of electrons from semileptonic decays reconstructed in the MC 
sample satisfying the selection criteria detailed in Section |5l and Nl^_^^ is the number 
of electrons from semileptonic decays produced in the signal process that satisfy the 
kinematic requirements of the cross-section definition using the MC information at the 
generator level. The kinematic variables and x at the true level were calculated using 
the four- momentum of the exchanged photon after possible initial-state radiation (ISR). 

The cross sections were corrected to the QED Born level, calculated using a running cou- 
pling constant, Oem, such that they can be compared directly to the NLO QCD predictions 
by HVQDIS. The radiative corrections were obtained using the Rapgap Monte Carlo as 
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Cr = CTBorn/crad, wheie (Trad IS the cross sectioii with full QED corrections (as used in the 
standard MC samples) and o"Born was obtained with the QED corrections turned off. The 
corrections are typically Cr ~ 1.05 rising to Cr ~ 1.10 for the high region. 

9 Systematic uncertainties 

The systematic uncertainties were calculated by varying the analysis procedure and then 
repeating the fit to the likelihood distributions |12]. The variations were made in a 
range such that the MC continued to provide a reasonable description of the data for the 
relevant distributions. The systematic uncertainties were determined bin by bin, unless 
stated otherwise. The main contributions came from the following sources, where the 
numbers in parentheses correspond to the uncertainty on the total cross section: 

1. DIS selection - the preselection cuts on the scattered electron were varied in both data 
and MC. The only cuts that had a significant effect were the cut on the energy, which 
was varied between 9 < Ee' < 11 GeV, the cut on the inelasticity, which was varied 
between 0.04 < ?/jb < 0.06, and the energy window for E — pz, which was varied by 
±4GeV(l}:^%); 

2. trigger efficiency - the uncertainty on the trigger efficiency was evaluated by comparing 
events taken with independent triggers (+1.2%); 

3. dE/dx simulation - both the mean and the width of the dE/dx distribution were 
varied in the MC separately and simultaneously by the uncertainty estimated from 
the data [23]. These two variations were then combined, giving a conservative estimate 
of the uncertainty on the dE/dx test function (lo;4%); 

4. tracking efficiency - the track-finding inefficiency in the data with respect to the 
MC was estimated to be at most 2%. The overall uncertainty due to this tracking 
inefficiency was determined by randomly rejecting 2% of all tracks in the MC and 
repeating the secondary- vertex finding (—3.4%); 

5. decay- length smearing - the fraction of events in the MC where the decay-length 
smearing was applied was varied by ±2% and the additional terms for the smearing 
of the tails were switched off {1l2^%)] 

6. p^' shape correction - the correction applied to the MC was switched off and increased 
by an additional 50% (zI'a^o); 

7. electron background - the relative contributions of the different electron sources in the 
MC were changed by varying separately the contributions from photon conversions, 
Dalitz decays, semileptonic decays from charm and DIS electrons by ±25% (^24%); 
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8. charm-spectrum reweighting - the correction to the c-decay electron spectrum in the 
MC using the CLEO data was varied by ±50% {tli%); 

9. energy scale - the global energy scale was varied in the MC by =f2% (l;j;'Q%); 

10. jet energy scale - the calorimetric part of the transverse jet energy in MC was varied 
by ±3% 

11. MC model dependence - the reweighting correction was varied by a factor of two 

12. electron reconstruction efficiency - the electron reconstruction efficiency in MC was 
varied by ±0.05 in the FCAL and RCAL regions and by =f0.05 in the BCAL region 

A series of further checks were made. The fit range was varied to check possible deficits 
in the background description. Selection cuts such as the Z vertex position or preselec- 
tion cuts such as on the dE/dx test function were varied before repeating the analysis. 
Another important check was the charge dependence. Separate fits were made for elec- 
tron and positron candidates for each lepton-beam charge separately as well as for the 
combined sample. All variations were found to be small and consistent with the expected 
fiuctuations due to statistics and were therefore not included in the systematic error. 

The individual contributions to the systematic uncertainties were added in quadrature, 
separately for the negative and the positive variations, to determine the systematic un- 
certainty of 177% for the total cross section. A ±2.0% overall normalisation uncertainty 
associated with the luminosity measurement was included in the uncertainty on the total 
cross section. 



10 Results 

The visible cross section for electrons from direct and indirect 6-quark decays with 0.9 < 
< 8 GeV in the range \r]^\ < 1.5 was measured in DIS events with > 10 GeV^ and 
0.05 < y < 0.7 and found to be 

CTfe^e = (71.8 ± 5.5(stat.)+5:5(syst.)) pb. 

This cross section includes all electrons and positrons from both b and b and no jet 
requirement was applied at the true level. This result can be compared to the HVQDIS 
NLO QCD prediction of 

= (6711?) pb. 
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where the uncertainty is calculated as described in Section HI This value agrees well with 
the measured cross section, which is a factor 1.3 higher than the Rapgap leading-order 
predictioijfl of 54.4 pb. This factor is used to scale the Rapgap predictions in Figs. [5] 
and [HI 

Differential cross sections as a function of and ?7e, and x are shown in Fig. 
Figure [6] shows the differential cross sections as a function of x, split into four different 
ranges. The figures also show the NLO QCD and the scaled Rapgap predictions. The 
cross-section values are given in Tables [IH31 Both the predictions from the NLO QCD 
calculations as well as the scaled Rapgap cross sections describe the data well. 



11 Extraction of F^^ 

The structure function Fg^^ can be defined in terms of the inclusive double-differential 
cross section (defined in analogy to Eq. [1]) as a function of x and Q^, 



dxdg2 xQ'^ 

where = 1 + {1 — yY and is the beauty contribution to the structure function F^. 

The electron cross section, cxfe^e, measured in bins of x and Q^, was used to extract 
at a reference point in the x-Q^ plane using 



dxdg2 dVfc^if /dxdQ2 ' 

where and dV^^if /dxdQ^ were calculated in the FFNS using the HVQDIS pro- 

gram. The uncertainty on the extrapolation from the measured range to the full kinematic 
phase space was estimated by varying the settings of the calculation (see Section |4]) for 
f!^'^^^ /{d^af}^^ /dxAQ"^) and adding the resulting uncertainties in quadrature. For each 
bin, a reference point in x and was defined (see Table H]) to calculate the structure 
function. The small correction for Ff^ is taken into account in the HVQDIS prediction. 

The structure function is shown in Fig. [7] as a function of x for nine different values 
of Q^. The values and the corresponding uncertainties are given in Table HJ To compare 
the result with previous measurements [31I11[T2], the earlier results were extrapolated to 
the values chosen in this analysis. For Q"^ > 10 GeV^, this measurement represents 
the most precise determination of F2'' by the ZEUS Collaboration. It is in good agree- 
ment with previous ZEUS analyses and the HI measurement. The NLO QCD prediction 



Note that the Rapgap predictions do not include the rcweighting correction discussed in Section|31 
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describes the data well. The same measurements are also shown as a function of for 
fixed X in Fig. [HI compared to several NLO and NNLO QCD predictions based on the 
fixed- or variable-flavour-number schemes [52]. For the HVQDIS prediction shown in this 
figure, the scale parametrisation /i = \\/ Q'^ + p\ + ml [3], was used. All the theoretical 
predictions shown provide a good description of the data. 

12 Conclusions 

Beauty production has been measured in DIS using semileptonic decays into electrons. A 
likelihood-ratio test function, adapted from a previous measurement, was used to identify 
the signal. The analysis benefited from the improved tracking in the HERA II data-set 
through the use of the measured decay length of weakly decaying B hadrons. 

The total cross section and differential cross sections as a function of x, Q^, and 77^ 
were determined. NLO QCD predictions calculated using the HVQDIS program describe 
the data well. The Rapgap Monte Carlo provides a good description of the shape of the 
differential distributions. 

The structure function F^^ was extracted from the double-differential cross section as 
a function of x and Q^. The measurement is in agreement with the results obtained 
from previous analyses using different techniques. For > 10 GeV^, this measurement 
represents the most precise determination of by the ZEUS Collaboration. The results 
were also compared to several NLO and NNLO QCD calculations, which provide a good 
description of the data. 
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Q2 
( GeV2) 


(pb/GeV2) (pb/GeV^) 


10 : 20 
20 : 40 
40 : 80 
80 : 200 
200 : 1000 


1.73±0.40l°:i 1.93l°i? 

1.05±0.18l°:J? 0.8410:13 
0.428±0.063l^:[}p 0.327lg:{]f? 
0.070±0.015l°:°?^ 0.087;[j:0}i 
0.0057±0.0014l0:0°?3 0.00661H00? 




X 


dab^Jdx da^i^^'/dx 
(pb) (pb) 


0.0002 : 0.0010 
0.0010 : 0.0020 
0.0020 : 0.0040 
0.0040 : 0.0100 
0.0100 : 0.1000 


34800±5700l^^|J|] 29700l^t|][] 
19400±2700l}^o° 14700;i[][] 
5800±1100l^[j|^ 5900l??[]o 
1200±310lio 1560;i[] 
38.4±12.ll|^ 48.5l|2 



Table 1: Differential cross sections for electrons from b-quark decays as a function 
of and x. The cross sections are given for > 10 GeV"^, 0.05 < y < 0.7, 
0.9 < < 8 GeV and jr^^l < 1.5. The first uncertainty is statistical and the second 
is systematic. In addition, the NLO QCD prediction and its uncertainty are given. 
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Pt 






(GeV) 


(pb/ GeV) 


(pb/ GeV) 


0.9 : 2.1 


36.9±6.ll^:? 


33.ll|i 


2.1 : 3.2 


12.2±2.0tJ:^ 




3.2 : 4.5 


3.08±0.90l°:^5 


A Qfi+0.59 
^•"^"-0.67 


4.5 : 8.0 


0.78±0.20l°:J^ 


n qc;+o.i3 



rf 




(pb) 


(pb) 


-1.5 : 


-0.5 


15.1±3.7l|^ 


13.411? 


-0.5 : 


0.0 


26.0±3.8l|^ 


26.711? 


0.0 : 


0.5 


30.3±5.ll^;^ 


30.0l^:e^ 


0.5 : 


1.5 


28.6±3.7lJ:^ 


oo 9+3.9 
^"-•■^-3.9 



Table 2: Differential cross sections for electrons from b-quark decays as a function 
of Pt ■ The cross sections are given for > 10 GeV'^, 0.05 < y < 0.7, 

0.9 < < 8 GeV and I'q'^l < 1.5. The first uncertainty is statistical and the second 
is systematic. In addition, the NLO QCD prediction and its uncertainty are given. 



( GeV^) 


dVb^e/dxdQ2 dVfe^if/dxdg^ 
(pb/GeV2) (pb/GeV2) 


10 : 20 0.0001 : 0.0004 
10 : 20 0.0004 : 0.0030 


2700±1200l?J][! 2500l|g[! 
300±100l^[! 480lJ°o 


20 : 60 0.0003 : 0.0012 
20 : 60 0.0012 : 0.0020 
20 : 60 0.0020 : 0.0060 


477±84l^^ 34311^ 
239±5ll^^ 180li° 
36±12l}^ 421^ 


60 : 400 0.0009 : 0.0035 
60 : 400 0.0035 : 0.0070 
60 : 400 0.0070 : 0.0400 


9.6±2.0lJ;^ 8.9lJ:° 
3.6±1.3lJ:° 5.01°;? 
0.23±0.12l°;?^ 0.47l°;gi 


400 : 1000 0.0050 : 0.1000 


0.013±0.010ino? 0.029in[]^ 



Table 3: Double- differential cross sections for electrons from b-quark decays as 
a function of x for four different ranges. The cross sections are given for 
> lOGeV^, 0.05 < y < 0.7, 0.9 < < 8 GeV and \r]^\ < 1.5. The first 
uncertainty is statistical and the second is systematic. In addition, the NLO QCD 
prediction and its uncertainty are given. 
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(GeV) X 




12 0.0002 
15 0.0013 
25 0.0005 
30 0.0013 
40 0.005 
80 0.002 
120 0.005 
180 0.013 
600 0.013 


0.0074 ±0.0033 l^Zr^SJ^ 

0021 +0 nnn? +0-0003 +o.ooo4 

U.UUZi ICU.UUU/ _o.0005 -0.0004 

1 52 +0 2 7 +0-0015 +0.0025 

U.UiOZ ICU.UUZ / _o.0019 -0.0029 

0110+0 0023 +0 0022 +0.0019 

U.UiiU ICU.UUZO _o.0017 -0.0021 

0.0041 ±0.0014 l°:ri°;r7 

0.0208 ±0.0043 

0110 +0 40 +0-0029 +0-0015 
U.UiiU inu.uu^u _o.ooi9 -0.0015 

0.0050 ±0.0027 

0.0089 ±0.0067 ^nS^^^S:^ 



Table 4: The structure function given for nine different values of Q"^ and x. 
The first error is statistical, the second systematic and the last is the extrapolation 
uncertainty. 
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Figure 1: Distribution of the likelihood-ratio test function for the electron hy- 
pothesis, T^E/dx' for e^^TT^^K'^^p and p. All histograms were normalised to unity. 
The selection cut at — 21nT^^yj^ < 3 is indicated by the vertical line. All other 
selection cuts were applied. 
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Figure 2: The distribution o/— 21nT, where T is the test function, using the 
beauty hypothesis for electron candidates, compared to the Monte Carlo expectation 
after the fit described in the text. The arrow indicates the region included in the fit 
(—2lnT < 10). The shaded areas show the fitted contributions for electrons from 
b-quark decays, electrons from other sources and the non-electron background. 
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Figure 3: Distributions of the variables related to the event selection, after applying 
all selection cuts: for the kinematic variables (a) Q"^ and (b) x, for (c) the energy 
of the scattered electron, Egi and (d) the transverse momentum of the electron 
candidate, p^. The variables E—pz and Z^tx, which were used for the event selection 
are shown in (e) and (f), respectively. The shaded areas show the MC expectations 
for the contributions for electrons from h- quark decays, electrons from other sources 
and the non- electron background as denoted in the figure, after applying the scale 
factors from the fit. The summed distribution is compared with the data distribution 
shown by the black points. 21 
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-2lnT<10 -2 In T< 1.5 




d/5d d/5d 

Figure 4: Distributions of (a) pJ^'-, (c) A(j) and (e) d/6d for all candidates that 
enter the fit satisfying — 21nT < 10. The same plots are shown for the beauty- 
enriched region (—2hiT < 1.5) in (b), (d) and (f). For details see the caption of 
Fig.m 
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Figure 5: Differential cross sections for electrons from b-quark decays as a function 
of the kinematic variables (a) Q"^ and (b) x, and the decay electron variables (c) 
and (d) rj^ . The cross sections are given for > 10 GeV^, 0.05 < y < 0.7, 
0.9 < < 8 GeV and \r]'^\ < 1.5. The measurements are shown as points. The 
inner error bar shows the statistical uncertainty and the outer error bar shows the 
statistical and systematic uncertainties added in quadrature. The solid line shows 
the NLO QCD prediction, with the uncertainties indicated by the band; the dashed 
line shows the scaled prediction from Rapgap. 
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Figure 6: Double- differential cross sections for electrons from b-quark decays as a 
function of x for different regions of . Other details as in the caption of Fig. 
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Figure 7: The structure function F^^ (filled symbols) as a function of x for 
nine different values of compared to previous results (open symbols). The inner 
error bars are the statistical uncertainty while the outer error bars represent the 
statistical, systematic and extrapolation uncertainties added in quadrature. The 
band represents the uncertainty on the NLO QCD prediction. Previous data have 
been corrected to the reference range of this analysis. 
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Figure 8: The structure function (filled symbols) as a function of for 
fixed values of x compared to previous results (open symbols). The inner error bars 
are the statistical uncertainty while the outer error bars represent the statistical, 
systematic and extrapolation uncertainties added in quadrature. The data have been 
corrected to the same reference x as the previous analysis f^jj. The measurements 
are compared to several NLO and NNLO QCD predictions (see text for details). 
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